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Using a data sample of 62 million Y(AS) — * BB decays collected between 1999 and 2001 by the 
BABAR detector at the PEP-II asymmetric-energy B Factory at SLAC we study events in which 
one neutral B meson is fully reconstructed in a final state containing a charmonium meson and 
the flavour of the other neutral B meson is determined from its decay products. The amplitude 
of the CP-violating asymmetry, which in the Standard Model is proportional to sin2/3, is 
derived from the decay time distributions. We measure sin2/3 = 0.75 ±0.09 (stat)±0.04 (syst) 
and |A| = 0.92 ± 0.06 (stat) ± 0.02 (syst). The latter is consistent with the Standard Model 
expectation of no direct CP violation. These results are preliminary. In addition, we report 
on precision measurements of the B lifetimes, and the B°-B° oscillation frequency Arrid- 



1 Introduction 

CP violation has been a central concern of particle physics since its discovery in 1964 in the decays 
of K° L decay!. An elegant explanation of the CP-violating effects in these decays is provided 
by the CP-violating phase of the three-generation Cabibbo-Kobayashi-Maskawa (CKM) quark- 
mixing matrix i. However, existing studies of CP violation in neutral kann decays and the 
resulting experimental constraints on the parameters of the CKM matrix eI do not provide a 
stringent test of whether the CKM phase describes CP violation i. In the CKM picture, large 
CP violating asymmetries are expected in the time distributions of B° decays to charmonium 
final states. 

In general, CP violating asymmetries are due to the interference between amplitudes with 
a weak phase difference. For example, a state initially produced as a B° (B°) can decay to 
a CP eigenstate such as J/ifiK® directly or can oscillate into a B° (B°) and then decay to 
J/ip Kg. With little theoretical uncertainty in the Standard Model, the phase difference between 



these amplitudes is equal to twice the angle (3 = arg [ —V^V^/V^V^ ] of the Unitarity Triangle. 
The CP-violating asymmetry in this mode allows a direct determination of sin2/3, and can thus 
provide a crucial test of the Standard Model. 

A B°B° pair produced in T(4S) decays evolves in a coherent P-wave until one of the B 
mesons decays. If one of the B mesons, referred to as B tag , can be ascertained to decay to a 
state of known flavour, i.e. B° or B°, at a certain time t tag , the other B, referred to as B rec , at 
that time must be of the opposite flavour as a consequence of Bose symmetry. Consequently, the 
oscillatory probabilities for observing B°B°, B°B° and B°B° pairs produced in T(4S) decays are 
a function of At = t rec — t tag , allowing mixing frequency and CP asymmetries to be determined 
if At is known. 

At the PEP-II asymmetric e + e~ collider H, resonant production of the T(4S) provides a 
copious source of B°B° pairs moving along the beam axis (z direction) with an average Lorentz 
boost of (/?7) = 0.55. Therefore, the proper decay-time difference At is, to an excellent ap- 
proximation, proportional to the distance Az between the two B°-decay vertices along the axis 
of the boost, At w Az/c(/3 r y). The average separation between the two B decay vertices is 
Az (Pj) ctb = 260 /im, while the RMS Az resolution of the detector is about 180 fim. 

The lifetime, mixing and sin2/3 analyses share a common analysis strategy: 

• select events where one B, labeled B rec , is fully reconstructed; 

• determine the vertex of the other B decay, B tSug , in the event by performing a vertex fit to 
the remaining charged particle trajectories, and compute At. 

At this point, one can perform an unbinned likelihood fit to the At distribution and determine 
the B lifetime. The mixing and CP asymmetry measurements require one additional step: 

• determine the flavour of the B tag decay. 

To determine the oscillation frequency Am^ events are selected where B rec is reconstructed in 
a neutral decay mode with a known flavour, such as D* + ir~ , and a simultaneous unbinned like- 
lihood fit to the At distributions of events where B Tec and Stag have opposite flavour (unmixed) 
and equal flavour (mixed) is performed. To determine CP asymmetries, B rec is a reconstructed 
CP final state such as J/ipK®, and the At distributions of events where B tag is a B° and B° 
respectively are fitted simultaneously. 

In order to establish the experimental technique, we first present precision measurements of 
the B° and B + lifetimes, and the B° B° oscillation frequency Am^. These measurements share 
the same vertexing algorithm and At determination as the measurement of the CP asymmetries, 
and, in case of the mixing measurement, the same flavour tagging algorithm. To eliminate 
possible experimenter's bias the parameter under study was hidden in all analyses until the 
event selection and reconstruction, fitting procedure and systematic errors were finalized. 

2 The BABAR detector and data sets 

The data used were recorded with the BABAR detector in the period October 1999-December 
2001. The total integrated luminosity of the data set is equivalent to 56 fb -1 collected near the 
T(4S) resonance. The corresponding number of produced BB pairs is estimated to be about 
62 million. The measurement of the charged and neutral B lifetimes is based on the initial 
23 million BB pairs, whereas the dataset used for the Am^ measurement includes the first 32 
million BB pairs. This latter sample was also used for the previously published measurement! 
of sin2/3; in contrast, the preliminary measurement of sin2/3 described here utilizes the entire 
data sample available. 

Since the BABAR detector is described in detail elsewhere I, only a brief description is given 
here. Surrounding the beam-pipe is a 5-layer silicon vertex tracker (SVT), which provides precise 



measurements of the trajectories of charged particles as they leave the e + e~ interaction point. 
Outside of the SVT, a 40-layer drift chamber (DCH) allows measurements of track momenta 
in a 1.5 T magnetic field as well as energy-loss measurements, which contribute to charged 
particle identification. Surrounding the DCH is a detector of internally reflected Cherenkov 
radiation (DIRC), which provides charged hadron identification. Outside of the DIRC is a 
CsI(Tl) electromagnetic calorimeter (EMC) that is used to detect photons, provide electron 
identification and reconstruct neutral hadrons. The EMC is surrounded by a superconducting 
coil, which creates the magnetic field for momentum measurements. Outside of the coil, the 
flux return is instrumented with resistive plate chambers interspersed with iron (IFR) for the 
identification of muons and long-lived neutral hadrons. We use the GEANT4 package i to 
simulate interactions of particles traversing the BABAR detector. 



3 Measurement of B Lifetimes and Am^ 

3.1 Exclusive B reconstruction 

The so-called B^ sample used for lifetime and mixing analyses consists of events where B rec 
is reconstructed in the modes B° -» D^~Tr + , £> w ~ p + , D^~af, J/i^K* and B + -» D^°tt + , 
J/t/jK + , t/j(2S)K + . Charged and neutral D* candidates are formed by combining a D° with a 
7T _ or 7T°. D° candidates are reconstructed in the decay channels K + n~ , K + tt~tt°, K + it + it~-k~ 
and KgTT + Ti~ and D~ candidates in the decay channels K + tt~tt~ and K®tt~. We reconstruct 
J/ip and ip(2S) in the decays to e + e~ and fj- + fi~ and the ip(2S) decay to J/ip ir + ir~ . 

Continuum e + e~ — ► qq background is suppressed by requirements on the normalized second 
Fox- Wolfram moment i for the event and on the angle between the thrust axes of B rec and 
of the other B in the event. B Tec candidates are identified by the difference AE between the 
reconstructed B energy and y/s/2 in the T(45) frame, and the beam-energy substituted mass 
tties calculated from yfs/2 and the reconstructed B momentum. We require ities > 5.2GeV/c 2 
and | AE| < 3o"ae- The distributions of mES for selected -Bfl av candidates is shown in Fig. |l|. 



Figure 1: Beam energy substituted mass distribution for selected B° (a) and B + (b) candidates. In 56 fb _1 , we 
reconstruct 15K B° and 13K B + flavour tagged signal events. Average signal purities for jties > 5.27GeV/c 2 , are 

85 % and 89 %, respectively. 
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3.2 At determination 



The decay time difference, At, between the two B decays is determined from the measured 
separation, Az = z rec — z tag , along the z axis between the reconstructed B rec (z rec ) and flavour- 
tagging decay -Bt ag vertex (zt ag ). This measured Az is converted into At with the known 
T(4S) boost, including a correction on an event-by-event basis for the direction of the B mesons 



with respect to the z direction in the T{AS) frame. The At resolution is limited by the z 
resolution of the tagging vertex. The -Btag decay vertex reconstruction starts from all tracks 
in the event except those incorporated in B rec . An additional constraint is provided by the 
calculated B tag production point and three-momentum, determined from the three-momentum 
of the fully reconstructed B TCC candidate, its decay vertex, and the average position of the 
interaction point (with a vertical size of 10 fim) and the T(4S) average boost. The derived B tag 
trajectory is fit to a common vertex with the remaining tracks in the event. Reconstructed Kg 
or A candidates are used as input to the fit in place of their daughters in order to reduce bias 
due to long-lived particles. Tracks with a large contribution to the x 2 are iteratively removed 
from the fit until all remaining tracks have a reasonable fit probability or all tracks are removed. 
For 99.5% of the reconstructed events the r.m.s Az resolution is 180 //m. 

Two different parameterizations are used to model the decay-time difference resolution. In 
the measurements of Am^ and sin2/?, the time resolution function is approximated by a sum of 
three Gaussian distributions (core, tail, and outlier) with different means and widths, 

(1) 

where S t = At — At true . For the core and tail Gaussians, the widths = S^x o At are the event- 
by-event measurement errors scaled by a common factor Sk ■ The scale factor of the tail Gaussian 
is fixed to the Monte Carlo value since it is strongly correlated with the other resolution function 
parameters. The third Gaussian, with a fixed width of 03 = 8ps, accounts for outlier events 
with incorrectly reconstructed vertices (less than 1% of events). The offsets hi are modeled to 
be proportional to a At, which is correlated with the weight that the remaining daughters of 
charm particles have in the tag vertex reconstruction. The tail and outlier fractions and the 
scale factors are assumed to be the same for all decay modes, since the precision of the -Btag 
vertex measurement is the limiting factor for the At resolution. This assumption is confirmed 
by Monte Carlo studies. 

The three Gaussian resolution function is less suited for the measurements of B lifetimes due 
to the large correlation between the resolution function parameters and the lifetimes, which leads 
to increased statistical errors. Studies with Monte Carlo simulations and data show that the sum 
of a zero-mean Gaussian distribution and its convolution with a one-sided exponential provides 
a good trade-off between statistical and systematic uncertainties in the lifetime measurement: 
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The parameters a are the fraction / in the core Gaussian component, a scale factor s for the 
per-event errors <taj, and the factor k in the effective time constant no At of the exponential 
which accounts for charm decays. At outlier events are modeled the same way as in the three 
Gaussian resolution function. The resolution functions differ only slightly between B° and B + 
mesons due to different mixtures of D~ and D° mesons in the Stag decays and we use a single 
set of resolution function parameters for both B° and B + in the lifetime fits. 



3.3 Lifetime results 

We extract the B + and B° lifetimes from an unbinned maximum likelihood fit to the At distri- 
butions of the selected B candidates. The probability for an event to be signal is estimated from 



toes fits (Fig. ||) and the toes value of the B Tec candidate. In the likelihood, the probability 
density for the signal events is given by 

/+oo 
e -|A*l/y(2 T )ft(A*- At',a At \a) d(At'), (3) 
-oo 

and the background At distribution for each B species is empirically modeled by the sum of a 
prompt component and a lifetime component convolved with the same resolution function, but 
with a separate set of parameters. The likelihood fit involves 17 free parameters in addition to 
the B° and the B + lifetimes: 12 to describe the background At distributions and 5 for the signal 
resolution function. The charged B lifetime t b + is replaced with t b + = r • t b o to estimate the 
statistical error on the ratio r = t b + / t b o . 

We determine the B and B + meson lifetimes and their ratio to be: 

t b o = 1.546 ± 0.032 (stat) ± 0.022 (syst) ps, 
t b + = 1.673 ± 0.032 (stat) ± 0.023 (syst) ps, and 
t b + /t b o = 1.082 ± 0.026 (stat) ± 0.012 (syst). 

These are the most precise published measurements to date0 and are consistent with the 
world averages 0. The resolution function parameters are consistent with those found in a 
Monte Carlo simulation that includes detector alignment effects. With the current data sample 
these measurements are still statistically limited. The dominant systematic errors arise from 
uncertainties in the description of the combinatorial background and of events with large At 
values, the use of a common time resolution function for B and B + and from limited Monte 
Carlo statistics. 

3.4 Flavour tagging 

After the daughter tracks of the B Tec are removed from the event, the remaining tracks are 
analyzed to determine the flavour of the -£>tag> and this ensemble is assigned a tag flavour, 
either B° or B°. For this purpose, flavour tagging information carried by primary leptons from 
semileptonic B decays, charged kaons, soft pions from D* decays, and more generally by high 
momentum charged particles is used to uniquely assign an event to a tagging category. 

Events are assigned a Lepton tag if they contain an identified lepton with a center-of- 
mass momentum greater than 1.0 or l.lGeV/c for electrons and muons, respectively. The 
momentum requirement selects mostly primary leptons by suppressing opposite-sign leptons 
from semileptonic charm decays. If the sum of charges of all identified kaons is non-zero, the 
event is assigned a Kaon tag. The final two tags involve a multi-variable analysis based on a 
neural network, which is trained to identify primary leptons, kaons, and soft pions, and the 
momentum and charge of the track with the maximum center-of-mass momentum. Depending 
on the output of the neural net, events are assigned either an NT1 (more certain) tag, an NT2 
(less certain) tag, or are considered not tagged (about 30% of events) and excluded from the 
analysis. The tagging power of the NT1 and NT2 tags comes primarily from slow pions, from 
kinematically recovering non-identified primary electrons and muons, and from kaons that do 
not pass the selection criteria for the Kaon category. 

Tagging assignments are made mutually exclusive by the hierarchical use of the tags. Events 
with a Lepton tag and no conflicting Kaon tag are assigned to the Lepton category. If no Lepton 
tag exists, but the event has a Kaon tag, it is assigned to the Kaon category. Otherwise the event 
is assigned to one of the two neural network categories. 

The effective tagging efficiency Qi = £j(l — 2wi) 2 , where £j is the fraction of events assigned 
to category i and Wi the probability of obtaining a wrong tag, is used as the basis for optimization 
of category selection criteria. The statistical errors on Ato^ and sin2/3 are proportional to 
where Q = J2Qi- The contributions of the various tagging categories to Q is shown in Table [j]. 



Table 1: Tagging efficiency e, average mistag fractions w, mistag differences Am = w(B°) — w(B°), and the 



derived Q (defined in the text) obtained from the likelihood fit to the -Bfl av and Bqp samples. 



Category 


£ 


(%) 


w 


(%) 


Aw ( c 


Yo) 


Q(%) 


Lepton 


11.1 


± 0.2 


8.6 


± 0.9 


0.6 ± 


1.5 


7.6 ± 0.4 


Kaon 


34.7 


± 0.4 


18.1 


± 0.7 


-0.9 ± 


1.1 


14.1 ± 0.6 


NT1 


7.7 


± 0.2 


22.0 


± 1.5 


1.4 ± 


2.3 


2.4 ± 0.3 


NT2 


14.0 


± 0.3 


37.3 


± 1.3 


-4.7 ± 


1.9 


0.9 ± 0.2 


All 


67.5 


± 0.5 










25.1 ± 0.8 



3.5 Mixing result 

The value of Ato^ is extracted from the tagged flavour-eigenstate B° sample with a simultaneous 
unbinned likelihood fit to the At distributions of both unmixed and mixed events. The PDFs 
for the unmixed (+) and mixed (— ) signal events for the i th tagging category are given by 

e -|At|/r 

H±{At,aAt\Arn d ,Wi,ai) = — [1 ± (1 - 2wi) cos Am d At] ®Tt(8 t ,cr/xt\ai). (4) 

4r 

Some resolution function parameters are allowed to differ for each tagging category to account 
for shifts due to inclusion of charm decay products in the tag vertex. The PDFs are extended 
to include background terms, different for each tagging category. The probability that a B° 
candidate is a signal event is determined from a fit to the observed toes distribution for its 
tagging category. The At distributions of the combinatorial background are described with 
a zero lifetime component and a non-oscillatory component with non-zero lifetime. Separate 
resolution function parameters are used for signal and background to minimize correlations. 

The At distributions of the signal (toes > 5.27GeV/c 2 ) overlaid with the projections of the 
likelihood fit, are shown in Fig. 0. In addition, the mixing asymmetry, 

(At) ^unmixedi. At) A^ m j :rea i(At) . . 

N un mixed( At) N m i xe( {(Al) 

is plotted. If flavour tagging and At determination were perfect, the asymmetry as a function 
of At would be a cosine with unit amplitude. 

The probability to obtain a likelihood smaller than that observed is 44%, evaluated with a 
parameterized Monte Carlo technique. The value of Ato^ obtained is 

Am d = 0.516 ± 0.016(stat) ± O.OlO(syst) ps -1 . (6) 

Since the parameters of the At resolution and the mistag rates w are free parameters in the fit, 
their contribution to the uncertainty on Am^ is included as part of the statistical error. The main 
contributions to the systematic errors are the choice of the signal At resolution description, its 
capability to handle outliers and various worst-case SVT misalignment scenarios (±0.005 ps -1 ), 
and by correlations between mistag rates and At resolution which are not explicitly modeled 
by the likelihood fit (±0.005 ps -1 ). Finally, the variation of the fixed B° lifetime within known 
errors leads to a systematic uncertainty of ±0.006 ps -1 . 

This is one of the single most precise mixing measurements availablfl, and is consistent 
with the current world average 0. 



4 Determination of sin2/3 



For the measurement of sin2/3, B rec is fully reconstructed in a CP eigenstate with eigenvalue 
VCP = -1 (J/ipK®, 4>(2S)K% or XaKf) or ±1 (J/ipK°), while B tSLg is tagged just as for the 
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Figure 2: Distributions of At for (a) unmixed and (b) mixed events in the signal region mes > 5.27GeV/c 2 . 
The data points are overlaid with the result from the fit, projected using the individual signal probabilities and 
event-by-event At resolutions, along with the simultaneously determined background distribution. Also shown 
(c) is the time-dependent mixing asymmetry „4 m i x (|At|) defined in the text. 

mixing measurement. The sample is further enlarged by including the mode J/t/)K*° (K*° — > 
Kgir ). However, due to the presence of even (L = 0, 2) and odd (L = 1) orbital angular 
momenta in the J/ipK*° system, there are rjcp = —1 and +1 contributions to its decay rate. 
These contributions are disentangled by incorporating their dependence on the transversity 
angles in each event into the likelihood frftl. The tties distributions (AE for J/ip K®) of the 
selected sample are shown in Fig. ||, and the detailed breakdown in Table ||. 

Table 2: Number of tagged events, signal purity, and result of fitting for CP asymmetries in the full CP sample 
and in various subsamples, as well as in the -Bfl av and charged B control samples. Purity is the fitted number of 



signal events divided by the total number of events 


in the AE and 


tties signal 


region defined in the text. 


Sample 


A^tag 


Purity 


(%) sin2/? 


Full CP sample 


1850 


79 


0.75 ± 0.09 


J/^K" S (Kg -> 7T+7T-) 


693 


96 


0.79 ± 0.11 


J/4jK° s (K° s ^TT°7T ) 


123 


89 


0.42 ± 0.33 


119 


89 


0.84 ± 0.32 


XciK° 3 


60 


94 


0.84 ± 0.49 


J/if>K° 


742 


57 


0.73 ± 0.19 


J/^K*° (K*° -> if°7r°) 


113 


83 


0.62 ± 0.56 


J/il>K», 1>(2S)K», XciK" only ( Vf - -1) 


995 


94 


0.76 ± 0.10 


Lepton tags 


176 


97 


0.73 ± 0.16 


Kaon tags 


504 


95 


0.75 ± 0.14 


NT1 tags 


117 


95 


0.86 ± 0.33 


NT2 tags 


198 


94 


0.84 ± 0.61 


B° tags 


471 


94 


0.79 ± 0.14 


B° tags 


524 


95 


0.73 ± 0.14 


-Bflav sample 


17546 


85 


0.00 ± 0.03 


Charged B sample 


14768 


89 


-0.02 ± 0.03 



The decay-time distribution of B decays to a CP eigenstate with a B° or B° tag can be 
expressed in terms of a complex parameter A that depends on both the B°-B° oscillation am- 
plitude and the amplitudes describing B° and B° decays to this final state Ej. The decay rate 
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Figure 3: Distribution of J71es for flavour tagged Bcp candidates selected in the final states a) J/ipKg (Kg — > 
7T+7T-), b) J/il>K% (K s -» ttV), c) V(2S)#<j, d) Xd^s, e) J/tpK*° {K*° -> KgTr ), and f) distribution of A£ 

for flavour tagged J/ip K° L candidates. 



JF+(.F_) when the tagging meson is a B°(B°) is given by 
^(AijCTAtl sm2(3,Wi,a,i) 

e -\At\/r B o 



1 ± (1 - 2wi 



21mX 



sin (Am^At) 



1 



1 + 



4r B o 

® K(Su<7At\ai), (7) 

where the + or — sign indicates whether the B tSLg is tagged as a B° or a B°, respectively. 

The distributions are much simpler when |A| = 1, which is the expectation of the Standard 
Model for decays like B° — > J/tp where all amplitudes which contribute to the decay have 
the same weak phase. In this particular case one is left with the phase difference introduced by 
mixing, i.e. A = rjcpe " . where rjcp is the CP eigenvalue of the final state. 

It is possible to construct a CP-violating observable which, neglecting resolution effects, is 
proportional to sin2/3: 

F+(At) -F-{At) 



Acp(At) 



.F + (At)+.F_(At) 



oc — rjcp(l ~ 2w) sin 2(3 sin Am^At. 
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Figure 4: Number of rjf = —1 candidates (J/ipKg, ij){2S)K®, XaK%) in the signal region a) with a B° tag N B o 
and b) with a B° tag ATg ; and c) the raw asymmetry (N B o — Ng )/(N B a + Ng ), as functions of At. The 
solid curves represent the result of the combined fit to the full Bcp sample. The shaded regions represent the 
background contributions. Figures d) - f) contain the corresponding information for the r)j = +1 mode J/ipK^. 
The likelihood is normalized to the total number of B° and B° tags. The value of sin2/3 is independent of the 
individual normalizations and therefore of the difference between the number of B° and B° tags. This difference 
is responsible for the small vertical shift between the data points and the solid curves. 



Since no time-integrated CP asymmetry effect is expected, an analysis of the time-dependent 
asymmetry is necessary. The interference between the two amplitudes, and hence the CP 
asymmetry, is maximal after approximately 2.1 B° proper lifetimes, when the mixing asym- 
metry goes through zero. However, the maximum sensitivity to sin2/3, which is proportional to 



-TAt 



sin Am^At, occurs in the region of 1.4 lifetimes. 



The value of sin2/3 can be extracted by maximizing the likelihood function 
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In C C p X 



ln:F + (Ai,<7At;sin2/3, Wi,a{,) + ^ ln.F_(Ai, cta*; sin2/3, Wi, cti, 



E 

B° tag B° tag 

(9) 

where the outer summation is over tagging categories i and the inner summations are over the 
B° and B° tags within a given tagging category. In practice, the fit for sin2/3 is performed 
on the combined i?fl av an d B>cp samples with a likelihood constructed from the sum of Eq. ||| 
and in order to determine sin2/3, the mistag fraction u>i for each tagging category, and the 
resolution parameters simultaneously. Additional terms are included in the likelihood to 
account for backgrounds and their time dependence. The determination of the mistag fractions 
and At resolution function for the signal is dominated by the high-statistics i?fl av sample. We 
fix t b o = 1.548 ps and Am^ = 0.472 ps -1 0. The largest correlation between sin2/3 and any 
linear combination of other free parameters is only 0.14. 



The simultaneous fit to all CP decay modes and the flavour decay modes yields: 



sin2/3 = 0.75 ± 0.09(stat) ± 0.04(syst) . (10) 

The dominant sources of systematic uncertainties are the choice of parameterization of the At 
resolution function, possible differences in the mistag fractions between the CP sample and 
the flavour sample, and uncertainties in the level, composition and CP asymmetry of the 
background in the selected events. The large sample of fully reconstructed events allows a 
number of consistency checks, including separation of the data by decay mode, tagging category 
and B tag flavour 

This analysis^ improves upon and supercedes the previously published result! It provides 
the single most precise measurement of sin2/? currently available and is consistent with the range 
implied by indirect measurements and theoretical estimates of the magnitudes of CKM matrix 
elements in the context of the Standard ModeS. 
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